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ABSTRACT 


Increasing nse of curved sandwich panels as aerospace structure components makes it vital to fully 
understand their thermostructural behavior and identify key factors affecting the open-mode debonding 
failure. Open-mode debonding analysis is performed on a family of curved honeycomb-core sandwich 
panels with different radii of curvature. The curved sandwich panels are either simply supported or 
clamped, and are subjected to uniform heating on the convex side and uniform cryogenic cooling on the 
concave side. The finite-element method was used to study the effects of panel curvature and boundary 
condition on the open-mode stress (radial tensile stress) and displacement fields in the curved sandwich 
panels. The critical stress point, where potential debonding failure could initiate, was found to be at the 
midspan (or outer span) of the inner bonding interface between the sandwich core and face sheet on the 
concave side, depending on the boundary condition and panel curvature. Open-mode stress increases with 
increasing panel curvature, reaching a maximum value at certain high curvature, and then decreases 
slightly as the panel curvature continues to increase and approach that of quarter circle. Changing the 
boundary condition from simply supported to clamped reduces the magnitudes of open-mode stresses and 
the associated sandwich core depth stretching. 

NOMENCLATURE 
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cross-sectional area of honeycomb cell (normal to honeycomb cell generatrix), in 2 

cross-sectional area of honeycomb cell wall (normal to the honeycomb cell 
generatrix), in 2 

Young’s modulus of face-sheet material, lb/in 2 

effective moduli of elasticity of honeycomb core, lb/in 2 

beam element for which the intrinsic stiffness matrix is given 

quadrilateral membrane element 

quadrilateral combined membrane and bending element 

shear modulus of face- sheet material, lb/in 2 

effective shear moduli of honeycomb core, lb/in 2 

depth of honeycomb core, in. 

joint location (grid point) of finite-element model 

half-span curved length along center line of curved sandwich panel, in. 

= 57.3 (/ /0), radius of curved sandwich panel, in. 

Structural Performance And Resizing finite-element computer program 
thermal protection system 

temperature of concave (lower) surface of curved sandwich panel, °F 
temperature of convex (upper) surface of curved sandwich panel, °F 
thickness of face sheets, in. 

coefficients of thermal expansion of face sheets, in/in-°F 
coefficients of thermal expansion of honeycomb core, in/in-°F 
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radial displacement of middle surface at midspan of curved sandwich panel, in. 

change of honeycomb core depth h (positive for stretching), in. 

maximum value of Ah, in. 

curvature angle (or half-span angle), deg 

Poisson’s ratio 

weight density of material used for face sheets and honeycomb core, lb/in 3 

effective weight density of honeycomb core structure, lb/m 3 

radial tensile stress (open-mode stress) in honeycomb core cell wall, lb/in 2 

tangential stress in concave side (lower) face sheet, lb/in 2 

tangential stress in convex side (upper) face sheet, lb/in 2 

effective radial tensile stress in honeycomb core, lb/in 2 

local maximum value of o r in a distribution of o r along core depth at a tangential 
station, lb/in 2 

maximum value of (o r ) d , lb/in 2 


Subscripts 

1, 2, 3 radial, tangential, and axial (longitudinal) directions 


INTRODUCTION 

Since the well-known successful aerospace application of sandwich structures in the all-wood- 
constructed British Mosquito fighter-bomber aircraft during World War II (refs. 1-4), the use of 
sandwich structural technology has become widespread in various aerospace structural applications (e.g., 
wings, tails, wall panels, webs of beams). The typical sandwich structure in panel form is fabricated with 
two relatively thin, high-strength face sheets separated by and bonded to opposite sides of relatively 
thick, low-density, low-strength core. The resulting sandwich structures are lightweight and have high 
flexural stiffness. 

In most aerospace applications, the sandwich panels are curved (e.g., fuselage glove of space shuttle 
orbiter, certain landing gear doors) or flat with constant core depth or variable core depths to form 
aerodynamic shapes (e.g., rotary wing blades, T-38 horizontal stabilizers). When applied to hypersonic 
flight vehicles such as space shuttle orbiter structures that are subjected to severe aerodynamic heating, 
the sandwich structures are protected with a thermal protection system (TPS) so that the structures can 
operate in wann temperatures with low thennal gradient across the sandwich core depth (ref. 5). If the 
thermal gradients across the core depth are too severe, the induced thermal moment could become strong 
enough to bend the sandwich panel and disturb the original panel shapes. Ko (ref. 6) studied this problem 
in great detail. 

The most extensively used sandwich structure in aerospace technology is the honeycomb-core 
sandwich structure. Because the honeycomb cell generatrix (a line whose motion generates a honeycomb 
cell wall) is perpendicular to the face sheets, only line contact ensures proper bonding between the face 
sheets and the honeycomb core. Corrosion or excess open-mode deformation (moving apart of two face 
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sheets because of bending) can cause the line-contact bonding to lose its integrity — a major drawback of 
the honeycomb sandwich structure. For titanium alloys, the so-called superplastically-formed/diffusion- 
bonding fabrication process could be used to form sandwich panels with a variety of core geometries and 
achieve diffusion-bonded surface contact between the face sheets and sandwich core. 

If the sandwich panel is curved, the mechanical or thermal bending in the direction of reducing the 
panel curvature could induce tensile stress (called open-mode stress) in the sandwich core depth direction 
because the two face sheets tend to separate from each other, causing potential debonding between the 
face sheets and sandwich core. For the solid curved beams (monolithic or laminated composites) 
subjected to open-mode bending, the open-mode stress distribution in the curved beam depth direction is 
arch- shaped; that is, zero values occur at the inner and outer boundaries and the maximum value occurs 
inward of the middle surface (refs. 7 through 15). The location of the maximum open-mode stress point 
moves toward the middle surface as the curved beam depth decreases (refs. 7, 8, 15). For a curved 
(horse- shoe- shaped or semi-elliptic- shaped) sandwich beam subjected to open-mode mechanical bending, 
however, distribution of the open-mode stress along the core depth direction is almost linear, with the 
value reaching a maximum at the inner bonding interface (between the sandwich core and face sheets), 
and tapering down slightly toward the outer bonding interface (ref. 16). 

If the open-mode bending is too severe, this unfavorable location of the maximum open-mode stress 
point in the curve sandwich panel will be the debonding crack nucleation site for inducing the catastrophic 
debonding failure of the panels. One of the past catastrophic failures of the curved sandwich panels 
occurred in the honeycomb landing gear door panels of a certain rotary-wing aircraft. During the landing 
approach when the landing gear doors were wide open, a strong gust of wind induced excess open-mode 
bending, and caused a catastrophic debonding failure, resulting in the total loss of one of the landing gear 
doors. Recently, the curved sandwich panels have been designed to reinforce the nozzle ramps of certain 
nonconventional rocket engines. In such application, the curved sandwich panels will be subjected to 
open-mode thermocryogenic bending because the convex side of the curved sandwich panels will be 
exposed to higher temperatures, and the concave side to cryogenic temperatures. 

Increasing use of curved sandwich panels as aerospace structure components makes it vital to 
fully understand their thermo structural behavior and identify key factors affecting the open-mode 
debonding failure. 

This report concerns the finite-element, open-mode debonding analysis of curved honeycomb-core 
sandwich panels subjected to thermocryogenic bending. The results of the analysis show how the open- 
mode stress distributions and sandwich panel deformations vary with the panel curvature and boundary 
condition. The information in this report could serve as guidelines in the effective design of failure-free 
curved sandwich panels that must function under thermocryogenic environment. 

DESCRIPTION OF PROBLEM 

The following sections describe the geometry of the family of curved sandwich panels, boundary 
conditions, and thermocryogenic loading condition used in the analysis. 

Geometry 

Figure 1 shows a cross-section of the curved honeycomb-core sandwich panel with core depth h, 
identical face-sheet thickness t s , half-span curve length / (measured along the center line), radius 
of curvature r, and the curvature angle (or half-span angle) 0. Panels with different curvatures are 
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generated by simply changing 0 and keeping h, /, and t s constant. Dimension r automatically changes 
with 0 according to r = 57.3(//0) because / is constant. Table 1 lists the dimensions of the analyzed 
sandwich panels. 


Table 1. Geometry of curved sandwich panels. 


h =1.4 in. 

/ = 31 in. 
t s = 0.032 in. 

0 = 90° (r = 19.74 in.). 

75° (r = 23.68 in.). 

60° (r = 29.60 in.). 

45° (r = 39.47 in.), 

30° (r = 59.21 in.). 

15° (r = 118.41 in.), 

5° (r = 355.23 in.), 

0° (r = oo, flat) 


The lineal - dimensions in table 1 (0 excluded) are similar to those of the curved honeycomb-core 
sandwich panels proposed for reinforcing the nozzle ramps of certain nonconventional rocket engines for 
future space transportation systems. In the range of 15° > 0 > 0°, panel deformed shapes for additional 
curvature angles 0 = 2.5°, 6.25°, 7.5°, 10°, and 12.5° also were examined to find the critical curvature 
angle 0 at which the panel deformed modal shape changes. 

Boundary Conditions 

The edges of the curved sandwich panels are either simply supported (fig. 2(a)) or clamped (fig. 2(b)). 
The method of simply supporting a sandwich panel is slightly different from the conventional way of 
simply supporting a solid panel. At the sandwich panel edges (fig. 2(a)), transverse rigid bars are attached 
and pin-joined to the sandwich edges, with the bar midpoints pinned to fixed points lying in the middle 
surface of the sandwich panel. This method of edge support allows (1) panel edges to rotate freely with 
respect to the horizontal axes (i.e., edges of middle surface), (2) panel edge depth to remain straight after 
deformation, and (3) maximum transverse shear deformation to take place at the panel edges. For the 
clamped boundary condition (fig. 2(b)), the two face sheets and sandwich core at the panel edges 
are clamped. 


Thermocryogenic Loading 

The curved sandwich panel is subjected to thermocryogenic loading. The convex side (upper surface) 
of the sandwich panel is uniformly heated to temperature T u = 280 °F, and the concave side (lower 
surface) uniformly cooled to cryogenic temperature T } = -320 °F. This temperature range is typical for 
laboratory testing the structural integrity of curved sandwich panels designed to reinforce the nozzle 
ramps of certain nonconventional rocket engines. Such thermocryogenic loading certainly induces 
open-mode bending, and raises concern about potential open-mode debonding failure of the curved 
sandwich panels. 


FINITE-ELEMENT ANALYSIS 

The following sections describe the finite-element modeling of the curved sandwich panels, and the 
material properties used for the face sheets and honeycomb core elements. 


4 



Finite-Element Modeling 


The open-mode linear elastic debonding analysis used the Structural Performance And Resizing 
(SPAR) finite-element computer code (ref. 17). For simplicity, only a segment of each curved 
honeycomb-core sandwich panel was considered. Thus, the problem became two dimensional. Because 
of symmetry with respect to the y-axis (fig. 1), only a semi-span of the curved sandwich panel segment 
was modeled. Each face sheet was modeled with one layer of isotropic quadrilateral combined membrane 
and bending E43 elements, and the sandwich core with eight layers of anisotropic quadrilateral membrane 
E41 elements. 

To simulate the pinned rigid rod at the panel edge (fig. 2(a)), each node at the edge of the semispan 
model was pin-connected to the fixed supporting point lying in the middle surface using beam E22 
element (for which the intrinsic stiffness matrix is given). The stiffness of the E22 elements was made 
very large to simulate the rigid rods. Temperature distribution along the sandwich core depth was 
assumed to be linear. Figure 3 shows a semi-span, finite-element model generated for curved sandwich 
panel segment with curvature angle 0 = 90°. Finite-element models for different panel curvatures were 
generated from the 0 = 90° model by simply changing the curvature angle 0 and the associated radius of 
curvature r without disturbing other dimensions. Thus, the total number of joint locations (JLOC) and of 
finite elements remained unchanged. Table 2 lists the size of the finite-element model for any curvature 
angle 0. 


Table 2. Size of finite-element model for any curvature angle 0. 


JLOC 2211 

E41 elements 1600 

E43 elements 400 

E22 elements (simply supported case only) 10 


Material Properties 

The material properties of the face sheets and the honeycomb core used for the finite-element models 
are of certain age -hardened steel, and are given, respectively, in tables 3 and 4. 


Table 3. Material properties of face sheets (age-hardened steel). 



70 °F 

280 °F 

-320 °F 

E, lb/in 2 

29.1 x 10 6 

28.05 x 10 6 

30.05 x 10 6 

G, lb/in 2 

10.4 x 10 6 

10.04 x 10 6 

10.6 xlO 6 

V 

0.302 

0.31 

0.285 

a, in/in- °F 

9.17 x 10 6 

9.17 x 10 6 

9.17 x 10~ 6 * 

p, lb/in 3 

0.287 

0.287 

0.287 


"Actual data not available. 
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Table 4. Effective material properties of 
honeycomb core (age-hardened steel). 


E n = 0.1 x 10 4 lb/in 2 

E 22 = 0.1 x 10 3 lb/in 2 

E 33 = 0.1 x 10 3 lb/in 2 

G 12 = 0.75 x 10 5 lb/in 2 

G 23 = 0.1 x 10 3 lb/in 2 

G 31 = 0.4 x 10 5 lb/in 2 

a i = = a 3 = 8.1 x 10 6 in/in-°F 

p hc = 0.432 x 10 2 lb/in 3 


In table 4, the subscripts {1, 2, 3) denote the radial, tangential, and axial (longitudinal) directions. 
Data for the effective Poisson’s ratios for the honeycomb core were not available at the time of analysis. 
The effective Poisson’s ratios for the honeycomb core are usually negligibly small (on the order of 10 3 to 
10 7 ; ref. 2); therefore, they are assumed to be zero in the analysis. 

The finite-element analysis computes the effective radial tensile stress a e ^ for the honeycomb core. 
This stress is then converted to the actual radial tensile stress o r (defined as open-mode stress) induced in 
the honeycomb core cell wall according to the following relationship: 


°r ~ C eff J 


= O 


eff 


Phc 


( 1 ) 


where A c and A w are, respectively, the cross-sectional areas of the honeycomb cell and cell wall (normal 
to the honeycomb cell generatrix); p and p hc are, respectively, the weight densities of the honeycomb 
core material and core structure. 


RESULTS 

This section presents the results of finite-element, open-mode analysis of curved honeycomb-core 
sandwich panels subjected to thermocryogenic bending, including panel deformations, core depth 
changes, open-mode stresses, and face-sheet stresses. 

Panel Deformations 

Figures 4 through 11 show the deformed shapes of the curved sandwich panels with different 
curvatures supported under different boundary conditions. The undeformed shapes are superimposed on 
the deformed shapes to show the relative positions of both. In all plots, the panel displacements are 
magnified for the sake of visualization, hr the figures, values of radial displacement 8 0 at the midspan of 
the middle surface of each sandwich panel is indicated. The sign of 8 0 is defined as positive for upward 
(outward) displacement and negative for downward (inward) displacement. For 8 = 15° only (fig. 9(a)), in 
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addition to the value of midspan radial displacement 8 0 (no longer maximum), the maximum radial 
displacements at two outer-span points are shown. 


For the simply supported panels with curvature angles varying from 0 = 90° to 0 =15° (figs. 4(a) 
through 9(a)), the midspan regions of the sandwich panels move inwardly (S 0 = negative), and the outer 
span regions bulge out. At 0 = 45°, 30°, 15° (figs. 7(a) through 9(a)), the panel deforms into shallow 
M shapes. For curvature angles 0 = 5° and 0 = 0°, (figs. 10(a), 11(a)), the panel bows upward deforming 
into arch shapes without cave-in regions. The causes of M-shaped and arch- shaped deformations will be 
discussed later in the section, “Face Sheet Stresses.” For the simply supported curved panels, the core 
depth stretching Ah (i.e., pulling apart of two face sheets) is maximum at the midspan, and gradually 
tapers down in tangential direction, and becomes zero at the panel supported edges because of constraint. 

For the clamped panels, the midspan radial S 0 is positive at 0 = 90° (fig. 4(b)) and then becomes 
negative as the curvature angle decreases (figs. 5(b) through 11(b)). At a curvature angle of 0 = 5° 
(fig. 10(b)), the midspan region appeal's cave-in because of the magnified displacement plot. When the 
panel is flat, 0 = 0°, (fig. 11(b)), the midspan downward displacement becomes infinitesimal. For the 
clamped curved panels (figs. 4(b) through 10(b)), the core depth stretching appeal's almost uniform over 
the entire span (except for the supported panel edges). As will be shown shortly, the core depth stretching 
Ah becomes maximum at the midspan or at the outer spans depending on the value of curvature angle 0. 

Figure 12 shows radial displacements S 0 at the midspan of the middle surface of the curved sandwich 
panel plotted as functions of curvature angle 0 for the two different boundary conditions. The simply 
supported case induced markedly larger magnitude of midspan displacements b 0 (a maximum of 1,357 
times larger at 0 = 0°) than the clamped case. For the simply supported case, the downward displacement 
of the midspan (8 0 = negative) is maximum at 0 = 90°, decreases monotonically with decreasing 0, turns 
to upward displacement (8 0 = positive) at approximately 0 = 14°, and finally increases steeply as 0 
approaches zero (flat panel). 

For the clamped case, the midspan displacement 8 0 is slightly upward (8 0 = positive) at curvature 
angle 0 = 90°, and turns downward (S 0 = negative) at around 0 = 85°. The downward displacement 
continues to increase with decreasing 0, reaching a peak at about 0 =10°, and then decreases to a very 
small negative value at 0 = 0°. 


Core Depth Changes 

Figure 13 shows the maximum sandwich core depth changes (A h) max plotted as functions of the 
curvature angle 0 for the two cases of boundary conditions. For the simply supported case, (A h) max is 
always at the midspan of the curved sandwich panels (flat panel excluded). For the clamped case, 
however, (A h) max is at the midspan only for the panels with low-curvature angles 45° > 0 > 5°. For the 
high-curvature angles 90° > 0 > 45°, the location of (A h) max shifts to the outer spans (near the edges) of 
the panel. 

The simply supported boundary condition induces higher values of (A h) max (a maximum of 
22 percent higher at 0 = 45°) than the clamped case. As the curvature angle 0 decreases from 0 = 90°, 
values of (A h) max increase slightly and reach their peak (indicated on the figures with downward arrows) 
at 0 = 75° for both boundary conditions, and then gradually decrease to very small negative values 
(contraction) at 0 = 0°. When the panel becomes flat (0 = 0°), no curvature effects can induce core 
stretching, and the sandwich core contracts slightly, because the cryogenic contraction of the sandwich 


7 



core depth overcomes the effect of its thermal expansion. Table 5 lists the numerical values of (A h) max 
used in plotting figure 13. 


Table 5. Maximum sandwich core depth changes (A h) max for different curvature angles 0. 


0, deg 

90 

75 

60 

45 

30 

15 

5 

0 

(A h) ma x 5 

Simply 

supported 

0.08736 

0.08949 

0.08838 

0.08198 

0.06413 

0.03497 

0.00841 

-0.00028 

in. 

Clamped 

0.07203* 

0.07383* 

0.07267* 

0.06629* 

0.05320 

0.03053 

0.01042 

-0.00028 


^Located at outer span; all others are located at midspan. 


Figure 14 shows the distribution of the sandwich core depth changes Ah (positive for stretching) for 
the simply supported case in the l-Q space, where l is defined in figure 1 . The distribution surface of Ah 
forms a distorted dome shape, with the apex at the midspan of 0 = 75° panel. The peak value of Ah is 
indicated with a downward arrow. For any curved sandwich panel, Ah is maximum at the panel midspan 
and gradually approaches zero toward the panel edge. The curve for connecting the midspan values of 
Ah (= (A h) max ) for different 0 is indicated in figure 14. This curve is identical to the (A h) max curve for 
the simply supported case shown in figure 13. At 0 = 0° (flat panel), the tangential distribution of Ah 
becomes flat with slightly negative values (contraction) over the entire span except for panel edges 
reflecting the overbalanced cooling of sandwich core as explained earlier. 

Figure 15 shows the distribution of sandwich core depth changes Ah in the /-0 space for the clamped 
case. The distribution surface of Ah for the clamped case is cascade shaped with the apex (peak value of 
Ah = (A h) max ) at the outer span of 0 = 75° panel as indicated with a downward arrow. For high -curvature 
angles 90° > 0 > 45°, the value of Ah tapers up slightly (almost linearly) from its midspan minimum value 
toward the outer span, reaching its peak before suddenly dropping to zero at the panel edges. For 
low-curvature angles 45° > 0 > 5°, Ah is maximum at the midspan and decreases infinitesimally and 
practically linearly toward the outer span and then drops suddenly to zero at the panel edges because of 
edge constraints. The (A h) max curve shown in figure 13 for the clamped case is actually a composite 
curve consisting of a segment of the curve at the outer span (curve with arrow sign attached) over the 
curvature angle range 90° > 0 > 45° where Ah = (A h) max , and a segment of the curve at the midspan over 
the range 45° > 0 > 5° where Ah = (A h) max (fig. 15). When the panel is flat (0 = 0°), Ah becomes slightly 
negative (contraction) and is practically constant over the entire panel span (except for the panel edges) 
because of overbalanced cooling as mentioned earlier. 

Open-Mode Stresses 

Figures 16 through 21 show the tangential distributions of local maximum open-mode stresses (o,.) rf 
induced in the curved sandwich panels with different curvature angles 0. As will be seen later, (o,.) d is at 
the inner bonding interface between the inner face sheet and the sandwich core. 
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For the simply supported case, the maximum open-mode stress (G r ) d = (o f ) is always at the 
midspan of the curved sandwich panel. The values of (<5 r ) d decrease monotonically (convex downward) 
in the tangential direction from its midspan maximum values (G r ) , and down to zero at the panel edges. 

For the clamped cases shown in the figures with dashed curves, the tangential distribution of (o r ) , 
remains almost constant over the span up to the panel edges where (o,.) d drops rapidly to zero because of 
clamping. The max im um open-mode stress (o r ) is at the midspan of the panel of low-curvature 
angle; 45° > 0 > 5°; and its location shifts to outer spans of the panel with high curvature angle 
90° > 0 > 45°. 

Figure 22 shows the plots of maximum open-mode stresses (G r ) max as functions of the curvature 
angle 0 for the two cases of boundary conditions. For the simply supported case, (G r ) is at the 
midspan of the inner bonding interface. For the clamped case (similar to (A h) max in figure 13), however, 
(o r ) is at the midspan of the inner bonding interface for the curvature range 45° > 0 > 5°, or at the 
outer spans of the inner bonding interface for the curvature angle range 90° > 0 > 45°. The simply 
supported case induces higher values of (o r ) than the clamped case, with a maximum of 32 percent 
higher at 0 = 90°. Tire peak values of (o r ) occur at 0 = 75° for both boundary conditions. The shapes 
of the (o r ) curves and the locations of peak values directly reflect characteristics of the (A h) max 
curves shown in figure 13. Table 6 lists the numerical values of (G r ) used in plotting figure 22. 


Table 6. Maximum open-mode stress (G r ) for different curvature angles 0. 


0, deg 

90 

75 

60 

45 

30 

15 

5 

1) 

(o,.) 

v > 'max 

lb/in 2 

Simply 

supported 

4,434 

4,533 

4,480 

4,134 

3,330 

1,940 

611 

0 

Clamped 

3,550* 

3,619* 

3,544* 

3,220* 

2,574 

1,489 

530 

0 


^Located at outer span; all others are located at midspan. 


Figure 23 shows the distribution of the local maximum open-mode stresses (G r ) d in the /-0 space for 
the simply supported case. Tire distribution of (o r ) d , like that of Ah (fig. 14), also forms a distorted dome 
shape, with its apex (peak value of (G r ) d = (G r ) ) indicated with a downward arrow, at the midspan 
of the 0 = 75° panel. For any curvature angle 0, the values of (o r ) d reach their peak (o r ) d = (G r ) max at 
the midspan, and gradually taper down to zero toward the panel edge. The curve for connecting the 
midspan values (o,.) is indicated in the figure. This curve is identical to the (o,.) curve for the 
simply supported case shown earlier in figure 22. At 0 = 0° (flat), (<5 r ) d is zero over the span. 

Figure 24 shows the local maximum open-mode stresses (o r ) d plotted in the /-0 space for the 
clamped case. The surface of ( o r ) y distribution for the clamped case looks similar to a cascade, with its 
apex (peak value of (<5 r ) d = (°,-) war ) at th e outer span of 0 = 75° panel (indicated on figure with a down- 
ward arrow). For any curvature angle 0, the value of (G r ) d appears almost constant over the panel span, 
and then suddenly drops to zero at the panel edge. When the panel is flat (0 = 0°), (o r ) rf is zero 
everywhere. The segment of curve at the outer span (curve with vertical arrow sign attached) over the 
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curvature angle range 90° > 0 > 45° where (o r ) d = (a r ) and the segment of curve at the midspan over 
the range 45° > 0 > 5° where (<5 r ) d = (<5 r ) max fonn the composite (a r ) cui’ve for the clamped case 
shown earlier in figure 22. 

Figures 25 and 26, respectively, show the radial distributions of open-mode stress o r along the y-axis 
(i.e., along the midspan core depth direction) for the simply supported and the clamped boundary 
conditions. For both types of boundary conditions, the radial distribution of o r for any curvature angle 
0 (0 * 0°) is almost linear, and tapers down very little toward the outer bonding interface. For the flat 
panel (0 = 0°), a r is zero everywhere in the core. 

For the simply supported case (fig. 25), c r = (G r ) is at the inner bonding interface (between the 
inner face sheet and the sandwich core), and the cui’ve indicated in the figure is identical to the (o r ) 
cui’ve shown earlier in figure 22 for the simply supported case. 

For the clamped case (fig. 26), the segment of (G r ) cui’ve lying in the region 0 < 45° is the 
midspan (o r ) cui’ve segment shown earlier in figure 22 for the clamped case. 

Face-Sheet Stresses 

Figures 27 through 34 show the tangential distributions of tangential stresses { Gg , Og ) induced in the 
concave (lower) and convex (upper) face sheets for different curvature angles 0 for the simply supported 
case. The plus (+) and minus (-) signs alongside the { Og , Gg ) curves imply tension and compression, 
respectively. The concave face sheet is under tension (Gg = positive) because of restrained cryogenic 
contraction, and the convex face sheet (except for the panel edge regions, figs. 29 through 33) is under 
compression (Og = negative) because of restrained thermal expansion. For a given curvature angle 0, 
the magnitude of Og is always larger than that of Og because of panel curvature and the unequal thermo- 
cryogenic loading condition. As 0 decreases from 0 = 90°, the difference between the magnitudes of { o e , 
Og } grows larger (figs. 27 through 33), and then becomes practically zero at 0 = 0° (fig. 34). Also, as 0 
decreases, Og in the convex face sheet near the panel edges begins to grow from negative to positive 
(figs. 29 through 33) because of the load transfer from the concave face sheet through the edge rigid rods. 
These panel edge tensile zones in the convex face sheet continue to grow larger as 0 becomes smaller, and 
finally extend to the entire panel span at 0 = 0° (fig. 34). 

The tangential stresses of opposite signs induced in the two face sheets (figs. 27 through 33) generate 
thermal moments [= (Og - Og )/?/ 2] that tend to bend the curved panels downward, thus appearing to be 
caved in. For a given curvature angle 0, the cave-in thermal moments reach maximum intensity at the 
midspan and gradually taper down toward the panel edges. At 0 = 0° (flat panel, fig. 34), the cave-in 
thermal moments diminish because both { Og , Gg } have practically identical small positive values that 
are constant over the entire panel span. 

Figure 35 shows the midspan magnitudes of { Gg , Gg } plotted as functions of curvature angle 0 for 
the simply supported case. The stress magnitude (Og-Og)/ 2 of the cave-in thermal moments 
[= (Gg - Og )h/2\ is also plotted. The value of (Og - Og)/2 increases as 0 decreases from 0 = 90°, 
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reaching a maximum at 0 = 15° where the M-shaped deformation is the most pronounced (fig. 9(a)), and 
then decreases rapidly toward zero at 0 = 0°. At 0 = 3°, Og changes sign from negative to positive and 
then reaches a small positive value practically identical to the value of Og (fig. 34) as 0 approaches 0 = 0°. 

Figure 36 shows the causes of the shallow M-shaped deformation of a typical 0 = 15°, (fig. 9(a)) 
simply supported curved panel. Because of the cryogenic contraction of the concave face sheet and the 
thermal expansion of the convex face sheet, the simply supported panel edges are forced to rotate about 
the fixed support points, causing upward bowing of the panel. Simultaneously, the induced cave-in 
thermal moments (maximum at the midspan) try to bend the panel downward, resulting in the M-shaped 
deformation. At lower curvature angles (e.g., 0 = 5°), the upward bowing effect begins to overshadow the 
cave-in bending effect, causing the curved panel to deform into a shallow arch shape (fig. 10(a)). The 
transition from M-shaped to arch- shaped deformation occurs in the neighborhood of 0 = 7° for the 
present study. 

Figures 37 through 39 show the tangential distributions of tangential stresses { Og , Og } induced in the 
concave and convex face sheets for different curvature angles 0 for the clamped case. The distributions of 
{ Og , Og }[or the cave-in thermal moments = (Og - Og )h/2] are almost constant over the panel span for 
the larger curvature angles between 0 = 90° and 0 = 30° (fig. 37), and starts to taper down slightly toward 
the panel edges at smaller curvature angles 0 = 15° to 0 = 5° (figs. 37 and 38). The differences in the 
magnitudes of { Og , Og } are almost inconspicuous for 0 = 90° to 0 = 30° panels, and start to grow larger 
with decreasing curvature angle 0. At 0 = 0° (fig. 39), the magnitudes of both { , Og ) are constant 

everywhere over the panel span (i.e., constant cave-in thermal moments, fig. 1 1 (b)). 

Figure 40 shows the midspan magnitudes of { Og , ojj } plotted as functions of curvature angle 0 for 
the clamped case. The stress magnitude (Og - Og )/2 for the cave-in thermal moment [= (Og - Og )h/2] is 
also plotted. The value of ( Og - Og )/2 increases as 0 decreases from 0 = 90°, reaching a maximum at 
0 = 5° where the most pronounced cave-in deformation occurs, (fig. 10(b)), and then decreases slightly as 
0 approaches 0 = 0°. 

Figure 41 shows the causes of cave-in deformation of a typical clamped panel (0 = 5°). Because of the 
clamped edges without rotations, upward bowing effect is constrained, and the cave-in thermal moments 
(almost constant over the entire panel span) bend the panel downward into the cave-in shape. 

DISCUSSION 

The near-linear radial distribution of o ; . across the core depth is typical for the curved honeycomb- 
core sandwich panels, and was observed also in horseshoe- and elliptic-curved honeycomb-core sandwich 
bars subjected to open-mode mechanical bending (ref. 15). For a classical solid curved beam on the other 
hand, the radial distribution of o r at any tangential cross-section is arch shaped, with zero values at the 
inner and the outer boundaries, and the local maximum value at a point slightly inward of the middle 
surface of the curved beam (refs. 7 through 15). 

For the present cases of curved sandwich panels (0 ^ 0°) under thermocryogenic loading, the inner 
bonding interface at the midspan or at the outer spans is, therefore, the potential debonding failure 
initiation region. The deformed shapes of curved sandwich panels presented (figs. 4 through 11) show that 
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the midspan and outer spans of the curved sandwich panels are the critical high deflection points. Thus, in 
using the curved sandwich panels as reinforcing structures operating under thennocryogenic environment, 
those critical high deflection points must be properly constrained (in addition to the panel edges) to 
preserve the original shapes. Because the clamped boundary condition induces smaller panel deflections 
and lower levels of open-mode stress field, by supporting the curved panel edges as close as the theoretical 
clamped condition, concerns of excess panel deflections and of the open-mode debonding failure could 
be m in imi z ed. 


CONCLUSIONS 

Finite-element open-mode debonding analysis was performed on a family of curved honeycomb-core 
sandwich panels subjected to thennocryogenic bending. The effects of panel curvature and boundary 
condition on the open-mode stress distributions and the defonnation fields were studied in detail. The key 
findings may be summarized in the following: 

1 . The peak panel deflection occurs at the midspan of the curved sandwich panel under both simply 
supported and clamped boundary conditions (exception: simply supported curved panel with 
curvature angle 0 = 15°, for which the peak panel deflections occur at outer spans). 

2. Tire radial distribution of open-mode stress in the core of a curved sandwich panel is practically 
linear, with local maximum and minimum values located, respectively, at the inner and outer 
bonding interfaces. The minimum value is only slightly lower than the maximum value. 

3. For a simply supported curved sandwich panel, the maximum open-mode stress point (or 
maximum core stretching point — the potential debonding failure initiation point) is always at the 
midspan of the inner bonding interface. 

4. For the clamped case, the maximum open-mode stress point (or maximum core stretching point) 
is at the midspan of the inner bonding interface for low-curvature panels only, and shifts to the 
outer spans of the inner bonding interface for high-curvature panels. 

5. The magnitude of the maximum open- mode stress (or maximum sandwich core stretching) 
increases with increasing panel curvature, reaching maximum at curvatures angle 0 = 75°, and 
then decreases slightly as the panel turns to a quarter circle (0 = 90°) under both simply supported 
and clamped boundary conditions. 

6. hi general, clamping a curved sandwich panel induces smaller panel deflections, lower open-mode 
stresses, and less core depth stretching than simply supporting the panel. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, February 3, 1999 
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Figure 1. Curved honeycomb-core sandwich panel subjected to heating and cryogenic cooling 
opposite sides. 




Figure 2. Two types of edge support conditions for the sandwich panel. 



















990045 

(a) Simply supported. 



990046 

(b) Clamped. 

Figure 7. Deformed shapes of curved honeycomb-core sandwich panel; T u = 280 °F; T t = -320 °F; 
0 = 45°. 



(a) Simply supported. 



(b) Clamped. 

Figure 8. Deformed shapes of curved honeycomb-core sandwich panel; T u = 280 °F; T t = -320 °F; 
0 = 30°. 
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(a) Simply supported. 



(b) Clamped. 

Figure 9. Deformed shapes of curved honeycomb-core sandwich panel; T u 
0 = 15°. 


280 °F; T, = -320 °F; 
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(a) Simply supported. 



990052 


(b) Clamped. 

Figure 10. Deformed shapes of curved honeycomb-core sandwich panel; T u = 280 °F; T t = -320 °F; 
0 = 5°. 
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(a) Simply supported. 



Figure 11. Deformed shapes of curved honeycomb-core sandwich panel; T u = 280 °F; T t = -320 °F; 
0 = 0° (flat). 



Figure 12. Plots of radial displacements b a (positive upward) at midspan of curved sandwich panel as 
functions of curvature angle 9; T u = 280 °F; T t = -320 °F. 
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Figure 13. Plots of honeycomb core maximum depth changes (A h) max (at midspan or outer span, positive 
for stretching) of curved sandwich panel as functions of curvature angle 0; T u = 280 °F; T t = -320 °F. 



Panel edge 


990058 

Figure 14. Distributions of honeycomb core depth change Ah (positive for stretching) in the /-0 space; 
T u = 280 °F; T f = -320 °F; shnply supported. 
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Figure 15. Distributions of honeycomb core depth change Ah (positive for stretching) in the /-0 space; 
T u = 280 °F; T, = -320 °F; clamped. 



Figure 16. Tangential distributions of local maximum open-mode stress (o ; .) (/ in curved honeycomb 
core; T u = 280 °F; T, = -320 °F; 0 = 90°. 
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Figure 17. Tangential distributions of local maximum open-mode stress in curved honeycomb 

core; T u = 280 °F; T, = -320 °F; 0 = 75°. 
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Figure 18. Tangential distributions of local maximum open-mode stress (o ; .) (/ in curved honeycomb 
core; 7V = 280 °F; T, = -320 °F; 0 = 60°. 




Figure 19. Tangential distributions of local maximum open-mode stress {G,.) d in curved honeycomb 
core; T u = 280 °F; T, = -320 °F; 0 = 45°. 



T,, = 280 °F 



0 = 30 


Figure 20. Tangential distributions of local maximum open-mode stress (o r ) (/ in curved honeycomb 
core; T u = 280 °F; T, = -320 °F; 9 = 30°. 
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T u = 280 °F 



e = i5° 


Figure 21. Tangential distributions of local maximum open-mode stress ( o r ) t in curved honeycomb core; 
T u = 280 °F; T, = -320 °F; 0 = 15°. 



Figure 22. Plots of maximum open-mode stress (o r ) . (at midspan or outer span) as function of curva- 
ture angle 9; T u = 280 °F; T, = -320 °F. 
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Figure 23. Distribution of local maximum open-mode stress (o r ) rf in the /-0 space; T u = 280 °F; 
T i = -320 °F; simply supported. 




| Peak value 


Figure 24. Distribution of local maximum open-mode stress (o r ) rf in the /-0 space; T u = 280 °F; 
Tj = -320 °F; clamped. 
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Figure 25. Radial distribution of open- mode stress o ; . along y-axis as a function of curvature angle 0; 
T u = 280 °F; T i = -320 °F; simply supported. 




Figure 26. Radial distribution of open-mode stress o r along y-axis as a function of curvature angle 0; 
T u = 280 °F; T, = -320 °F; clamped. 


29 




990072 


Figure 27. Tangential distributions of tangential stresses {Oq, Oq} in the face sheets; T u = 280 °F; 
T i = -320 °F; 0 = 90°; simply supported. 
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Figure 28. Tangential distributions of tangential stresses {g’q, Oq} in the face sheets; T u = 280 °F; 
T j = -320 °F; 0 = 75°; simply supported. 
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Figure 29. Tangential distributions of tangential stresses {Og, Og} in the face sheets; T u = 280 °F; 
T i = -320 °F; 0 = 60°; simply supported. 
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Figure 30. Tangential distributions of tangential stresses {Og, Og} in the face sheets; T u = 280 °F; 
T j = -320 °F; 0 = 45°; simply supported. 


31 




990075 


Figure 31. Tangential distributions of tangential stresses {Cq, Oq} in the face sheets; T u = 280 °F; 
T ! = -320 °F; 0 = 30°; simply supported. 



Figure 32. Tangential distributions of tangential stresses {Oq, Oq} in the face sheets; T u = 280 °F, 
T i = -320 °F; 0 = 15°; shnply supported. 
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Midspan Spanwise station Panel edge 

990077 

Figure 33. Tangential distributions of tangential stresses {Oq, Oq} in the face sheets; T u = 280 °F; 
T i = -320 °F; 0 = 5°; simply supported. 
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Figure 34. Tangential distributions of tangential stresses {Oq, Oq} in the face sheets; T u = 280 °F; 
T j = -320 °F; 0 = 0°; simply supported. 
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Figure 35. Face sheets tangential stresses {Og, Og} at panel midspan as functions of curvature angle 0; 
T u = 280 °F; T i = -320 °F; simply supported. 
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Figure 36. M-shaped deformation of the simply supported curved sandwich panel caused by the 
combined effect of the cave-in thermal moments induced in the panel and the bowing rotations at the 
panel edges; T u = 280 °F; T, =-320 °F; 0 = 15°. 
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Midspan Spanwise station Panel edge 
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Figure 37. Tangential distributions of tangential stresses {Oq, Oq) in the face sheets; T u = 280 °F; 
T i = -320 °F; 9 = 90°~15°; clamped. 
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Figure 38. Tangential distributions of tangential stresses {Oq, Oq} in the face sheets; T u = 280 °F; 
Tj = -320 °F; 9 = 5°; clamped. 
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Figure 39. Tangential distributions of tangential stresses {cjj, Oq} in the face sheets; T u = 280 °F; 
T i = -320 °F; 0 = 0°; clamped. 
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Figure 40. Face sheets tangential stresses { c’q, Oq } at panel midspan as functions of curvature angle 0; 
T u = 280 °F; T f = -320 °F; clamped. 




36 




Figure 41. Cave-in deformation of the clamped curved sandwich panel caused by the cave-in thermal 
moments induced in the panel; T u = 280 °F; T t = -320 °F; 9 = 5°. 
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